Phospho-seq: Integrated, multi-modal profiling of intracellular protein dynamics in single cells.

Cell signaling plays a critical role in regulating cellular behavior and fate. While multimodal single-cell sequencing technologies are rapidly advancing, scalable and flexible profiling of cell signaling states alongside other molecular modalities remains challenging. Here we present Phospho-seq, an integrated approach that aims to quantify phosphorylated intracellular and intranuclear proteins, and to connect their activity with cis-regulatory elements and transcriptional targets. We utilize a simplified benchtop antibody conjugation method to create large custom antibody panels for simultaneous protein and scATAC-seq profiling on whole cells, and integrate this information with scRNA-seq datasets via bridge integration. We apply our workflow to cell lines, induced pluripotent stem cells, and 3-month-old brain organoids to demonstrate its broad applicability. We demonstrate that Phospho-seq can define cellular states and trajectories, reconstruct gene regulatory relationships, and characterize the causes and consequences of heterogeneous cell signaling in neurodevelopment.

CPA-Perturb-seq: Multiplexed single-cell characterization of alternative polyadenylation regulators.

Most mammalian genes have multiple polyA sites, representing a substantial source of transcript diversity that is governed by the cleavage and polyadenylation (CPA) regulatory machinery. To better understand how these proteins govern polyA site choice we introduce CPA-Perturb-seq, a multiplexed perturbation screen dataset of 42 known CPA regulators with a 3' scRNA-seq readout that enables transcriptome-wide inference of polyA site usage. We develop a statistical framework to specifically identify perturbation-dependent changes in intronic and tandem polyadenylation, and discover modules of co-regulated polyA sites exhibiting distinct functional properties. By training a multi-task deep neural network (APARENT-Perturb) on our dataset, we delineate a cis-regulatory code that predicts responsiveness to perturbation and reveals interactions between distinct regulatory complexes. Finally, we leverage our framework to re-analyze published scRNA-seq datasets, identifying new regulators that affect the relative abundance of alternatively polyadenylated transcripts, and characterizing extensive cellular heterogeneity in 3' UTR length amongst antibody-producing cells. Our work highlights the potential for multiplexed single-cell perturbation screens to further our understanding of post-transcriptional regulation in vitro and in vivo.

Genomics analysis of Drosophila sechellia response to Morinda citrifolia fruit diet.

Drosophila sechellia is an island endemic host specialist that has evolved to consume the toxic fruit of Morinda citrifolia, also known as noni fruit. Recent studies by our group and others have examined genome-wide gene expression responses of fruit flies to individual highly abundant compounds found in noni responsible for the fruit's unique chemistry and toxicity. In order to relate these reductionist experiments to the gene expression responses to feeding on noni fruit itself, we fed rotten noni fruit to adult female D. sechellia and performed RNA-sequencing. Combining the reductionist and more wholistic approaches, we have identified candidate genes that may contribute to each individual compound and those that play a more general role in response to the fruit as a whole. Using the compound specific and general responses, we used transcription factor prediction analyses to identify the regulatory networks and specific regulators involved in the responses to each compound and the fruit itself. The identified genes and regulators represent the possible genetic mechanisms and biochemical pathways that contribute to toxin resistance and noni specialization in D. sechellia.

The CAR-mRNA Interaction Surface Is a Zipper Extension of the Ribosome A Site.

The ribosome CAR interaction surface behaves as an extension of the decoding center A site and has H-bond interactions with the +1 codon, which is next in line to enter the A site. Through molecular dynamic simulations, we investigated the codon sequence specificity of this CAR-mRNA interaction and discovered a strong preference for GCN codons, suggesting that there may be a sequence-dependent layer of translational regulation dependent on the CAR interaction surface. Dissection of the CAR-mRNA interaction through nucleotide substitution experiments showed that the first nucleotide of the +1 codon dominates over the second nucleotide position, consistent with an energetically favorable zipper-like activity that emanates from the A site through the CAR-mRNA interface. Moreover, the CAR/+1 codon interaction is affected by the identity of nucleotide 3 of +1 GCN codons, which influences the stacking of G and C. Clustering analysis suggests that the A-site decoding center adopts different neighborhood substates that depend on the identity of the +1 codon.

Arginine Methylation Regulates Ribosome CAR Function.

The ribosome CAR interaction surface is hypothesized to provide a layer of translation regulation through hydrogen-bonding to the +1 mRNA codon that is next to enter the ribosome A site during translocation. The CAR surface consists of three residues, 16S/18S rRNA C1054, A1196 (E. coli 16S numbering), and R146 of yeast ribosomal protein Rps3. R146 can be methylated by the Sfm1 methyltransferase which is downregulated in stressed cells. Through molecular dynamics analysis, we show here that methylation of R146 compromises the integrity of CAR by reducing the cation-pi stacking of the R146 guanidinium group with A1196, leading to reduced CAR hydrogen-bonding with the +1 codon. We propose that ribosomes assembled under stressed conditions have unmethylated R146, resulting in elevated CAR/+1 codon interactions, which tunes translation levels in response to the altered cellular context.